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ABSTRACT 

A nonlinear two dimensional fluid model of whistler turbulence is developed that non- 
linearly couples wave magnetic field with electron density perturbations. This coupling 
leads essentially to finite compressibility effects in whistler turbulence model. Interest- 
ingly it is found from our simulations that despite strong compressibility effects, the 
density fluctuations couple only weakly to the wave magnetic field fluctuations. In a 
characteristic regime where large scale whistlers are predominant, the weakly coupled 
density fluctuations do not modify inertial range energy cascade processes. Conse- 
quently, the turbulent energy is dominated by the large scale (compared to electron 
inertial length) eddies and it follows a Kolmogorov-like k~ 7 ^ 3 spectrum, where k is a 
characteristic wavenumber. The weak coupling of the density fluctuations is explained 
on the basis of a whistler wave parameter that quantifies the contribution of density 
perturbations in the wave magnetic field. 

Key words: (magnetohydrodynamics) MHD, (Sun:) solar wind, Sun: magnetic fields, 
ISM: magnetic fields 



1 INTRODUCTION 

Whistler wave regime is ubiquitously present in many space, 
astrophysical and laboratory plasmas. For instance, in the 
magnetospheric plasma, electrons can be accelerated by 
whistler-mode and compressional ULF (fast mode waves) 
turbulences near the Earth's synchronous orbit (Li et al, 
2005). The whistler- mode turbulence can accelerate sub- 
storm injection electrons with several hundreds of keV 
through wave-particle gyroresonant interaction and hence 
may play an important role in the electron acceleration dur- 
ing substorms (Li et al, 2005). Vetoulis & Drake (1999) de- 
scribe whistler turbulence at the magnetopause. Whistler 
mode turbulence can be triggered by electron beams in 
earth's bow shock (Tokar et al, 1984). In the solar wind 
plasma, observations have identified a spectral break in the 
solar wind magnetic field spectrum (Goldstein et al 1995, 
Leamon et al 1999). The mechanism leading to the spectral 
break has been unclear and thought to be either mediated 
by the kinetic Alfven waves (Hasegawa 1976), or by electro- 
magnetic ion-cyclotron- Alfven waves (Wu & Yoon, 2007), 
or whistler cascade regime (Gary et al, 2008), or by a class 
of fluctuations that can be dealt within the framework of 
the Hall magnetohydrodynamic plasma model (Alexandrova 
et al 2007, 2008; Shaikh & Shukla 2008, 2009). Stawicki et 
al (2001) argue that Alfven fluctuations are suppressed by 
proton cyclotron damping at intermediate wavenumbers so 
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the observed power spectra are likely to consist of weakly 
damped magnetosonic and/or whistler waves which are dis- 
persive unlike Alfven waves. Moreover, turbulent fluctua- 
tions corresponding to the high frequency and kpi ^> 1 
regime (where k is wavenumber, and pi is ion gyroradius) 
lead to a decoupling of electron motion from that of ion such 
that the latter becomes unmagnetized and can be treated as 
an immobile neutralizing background fluid. While whistler 
waves typically survive in the higher frequency (and the 
corresponding smaller length scales) part of the solar wind 
plasma spectrum, their role in influencing the inertial range 
turbulent spectral cascades is still debated (Biskamp et al, 
1996; Shaikh & Zank, 2003, 2005; Shaikh & Shukla, 2009, 
2008). 

Biskamp et al (1999) performed two and three dimen- 
sional simulations of incompressible (density perturbations 
are ignored) electron MHD model to demonstrate that the 
energy spectrum follows a k~ 5 ^ 3 law for kd e > 1 and k~ 7 ^ 3 
for kd e < 1. They further reported that the 3D spec- 
tral properties are similar to those in 2D. This was lately 
confirmed by Shaikh (2009) using 3D simulations. Cho & 
Lazarian (2004) performed 3D simulations of incompress- 
ible electron MHD model to study anisotropic scaling that 
relates the parallel and perpendicular wavenumbers through 
fen ~ k 1 / 3 ■ In a much detailed work, Cho & Lazarian (2009) 
examined the anisotropy in the electron MHD and showed 
that the high-order statistics in electron MHD admit a scal- 
ing that is similar to the She-Leveque scaling for incompress- 
ible hydrodynamic turbulence. 
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While there exists considerable literature describing 
the anisotropic cascades in whistler turbulence, the role of 
whistler wave in the spectral transfer of energy is still de- 
bated (Biskamp et al, 1996; Shaikh & Zank, 2003, 2005; 
Shaikh & Shukla, 2009, 2008). For instance, the Kolmogorov 
like dimensional arguments indicate that propagation of 
whistlers in the presence of a mean or an external constant 
magnetic field may change the spectral index of the inertial 
range turbulent fluctuations from k~ 7 ^ 3 to fc -2 (Biskamp et 
al, 1996). By contrast, the numerical simulations (Biskamp 
et al, 1996; Shaikh & Zank, 2003, 2005) suggest that whistler 
waves do not influence the spectral migration of turbulent 
energy in the inertial range despite strong wave activity 
and that the turbulent spectra corresponding to the electron 
fluid fluctuations in whistler turbulence continue to exhibit 
a Kolmogorov-likc fc~ 7 / 3 spectrum. 

What is not clear from these work is the quatitative 
role of whistler and the corresponding mode coupling in- 
teractions that mediate the inertial range turbulent spectra. 
Furthermore, much of the work described above (Biskamp et 
al, 1996, 1999; Shaikh & Zank, 2003, 2005; Shaikh & Shukla, 
2009, 2008; Cho & Lazarian, 2004, 2009) ignore the effect 
of density perturbations on the whistler mode turbulence. It 
is unclear if density fluctuations in the electron fluid mod- 
ify the turbulent cascade properties. Since density fluctua- 
tions are critically important in many space and laboratory 
plasma phenomena, their role in whistler turbulence needs 
to be investigated. 

The central object of this paper is to explore the non- 
linear turbulent processes mediated by whistler waves in the 
presence of density perturbations. We will investigate non- 
linear turbulent fluctuations, based on nonlinear fluid sim- 
ulations, in u> > uj ci regime where correlation length scales 
of turbulence are comparable to the electron inertial length 
scales. Understanding of whistler turbulence in the presence 
of density fluctuations is important in the context of so- 
lar wind plasma (Krafft & Volokitin, 2003; Saito et al, 2008; 
Stawicki et al, 2001; Gary et al, 2008; Ng et al, 2003; Vocks et 
al, 2005; Salem et al, 2007; Bhattacharjee et al, 1998), mag- 
netic reconnection in the Earth's magnetosphere (Wei et al, 
2007) to interstellar medium (Burman, 1975) and astrophys- 
ical plasmas (Roth, 2007) where characteristic fluctuations 
can typically be of several astronomical units. These are only 
a few of the numerous other studies. For more literature, 
the readers can refer to the simulation work by Biskamp 
et al (1996) and others including Shukla (1978), Shukla et 
al (2001), Galtier (2008), Urrutia et al (2008), Saito et al 
(2008), Bengt & Shukla (2008), Shaikh (2009), Cho & Lazar- 
ian (2004, 2009) and numerous references therein. 

This paper is organized as follows. Section 2 describes 
nonlinear whistler wave model for finite density perturba- 
tion. The finite density perturbations are described in terms 
of plasma beta (ratio of magnetic energy and electron fluid 
pressure). Section 3 deals with nonlinear fluid simulations. 
Energy spectra are discussed in Section 4 and section 5 de- 
scribes the effect of the density fluctuations in whistler tur- 
blence. In section 6, we discuss spectral anisotropy in in- 
homogeneous whistler turbulence. Section 7 concludes our 
work. 



2 WHISTLER WAVE MODEL 

Fluctuations in a magnetized plasma excite whistler modes 
when they propagate along a mean or background mag- 
netic field with characteristic frequency uj > Lu ci and the 
length scales are c/oj pi < £ < c/uj pe , where uj pi ,ui ve are the 
plasma ion and electron frequencies. In such a high frequency 
regime, the ions do not have time to respond to the elec- 
tron motions. Hence the electron dynamics plays a critical 
role in determining the nonlinear interactions while the ions 
merely provide a stationary neutralizing background against 
fast moving electrons and behave as scattering centers. The 
whistler turbulence can be described by the electron mag- 
netohydrodynamics (EMHD) model of plasma (Kingsep et 
al, 1990) that deals with the single fluid description of quasi 
neutral plasma. The EMHD model has been discussed in 
considerable detail in earlier work (Kingsep et al, 1990; 
Biskamp et al, 1996, 1999; Dastgeer et al, 2000a,b; Shaikh & 
Zank, 2003, 2005; Cho & Lazarian, 2004, 2009). In whistler 
modes, the currents carried by the electron fluid are im- 
portant, and we therefore write down only those equations 
which are pertinent to electron motion. These are electron 
fluid momentum, electric field, currents, and electron conti- 
nuity equations, 
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Here V • P = VP + V • II, the sum of pressure and stress 
tensors. The pressure becomes highly anisotropic in presence 
of a strong background magnetic field, especially in the low 
beta solar wind plasma. Hence the total pressure consists of 
the isotropic (VP) and the anisotropic (V ■ II) parts. The 
remaining equations are B = V x A, J = — en e V e , V B = 0. 
Here m e ,n e , V e are the electron mass, density and fluid ve- 
locity respectively. E, B respectively represent electric and 
magnetic fields and <f>, A are electrostatic and electromag- 
netic potentials. The remaining variables and constants are, 
the collisional dissipation /i, the current due to electrons flow 
J, and the velocity of light c. The displacement current in 
Ampere's law Eq. (3) is ignored because u) ce /Lu pe < 1. The 
plasma is assumed to be quasineutral, hence n e w m = n. 
Density fluctuations are considered to be incompressible at 
the leading order, but a first order compressibility is included 
to describe finite electron plasma beta (beta is the ratio of 
plasma pressure and magnetic field energy) effects. We adopt 
the approach described by Abdalla et al (2001) to include 
the finite density perturbation effect in the whistler wave 
model that couples density with pressure perturbations. The 
density field is determined from Poisson's equation and mo- 
mentum equation as follows, 



-4nn e e = V ■ E : 



— V - (V e x B) 

c 



(-) 

\n e e/ 



(5) 



where electron inertia and pressure contribution is negligi- 
bly small. Here perturbed density (n e ) is small compared 
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to the mean background density (no) such that n e /no < 1. 
Furthermore, lengthscale on which density is varying (Xn) is 
small compared to the characteristic length (£) i.e £xn < 1). 
The leading order electron fluid velocity can then be associ- 
ated with the rotational magnetic field through 



-V x B. 



(6) 



4irnoe 

Substitution of Eq. (6) into Eq. (5) with further simplifica- 
tion gives 



+V 2 /3 

Bo 



(7) 



where z is the direction of the background magnetic field 
and ft = 4-kp/Bq. The above relation is further consistent 
with n e /no < 1 since uj ce /uj pe < 1 and B/B < 1 in the low 
beta whistler mode turbulence. This inequality may how- 
ever not hold for the kd e S> 1 modes. In any event, the 
latter is inaccessible by fluid theory as kinetic effects begin 
to play a crucial role for the smallest scales. We nonetheless 
restrict ourselves to kde > 1 only where the characteristic 
length scales are marginally smaller than electron inertial 
length scales. On taking the curl of Eq. (1) and, after slight 
rearrangement of the terms, we obtain 



dSl' 
~dt 

where 



+ V X (V e X 0') = V X ( -— ) -flVxV 

V mn e J 



(8) 



n' = v x p = d 2 e v 2 B - b. 

It can be seen from Eq. (8) that in the ideal whistler mode 
turbulence (i.e. neglecting the term associated with the 
damping /i), the Curl of generalized electron momenta is 
frozen in the electron fluid velocity. This feature is strik- 
ingly similar to Alfvenic turbulence where the magnetic field 
is frozen in the ideal two fluid plasma (Biskamp, 2003). Us- 
ing electron continuity equation Eq. (4) in combination with 
Eq. (8), we obtain 



u)t 



+ V e - V 



n' 



— • V V e + V X 
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+uV x V e 



(9) 



We next introduce normalized generalized vorticity ft as fol- 
lows. 

n e u) ce B Bo \n / -Bo 

On substituting fl' into Eq. (8) and using appropriate vec- 
tor identities, we obtain the three-dimensional normalized 
equation of EMHD describing the evolution of the magnetic 
field fluctuations in whistler wave, 



V x ^1 +M?V 2 B. (10) 
V ran 



^ + V e ■ Vf2 - n • VV e 

at 

The length scales in Eq. (10) are normalized by the electron 
skin depth d e — c/uj pe i.e. the electron inertial length scale, 
the magnetic field by a typical amplitude Bo, and time by 
the corresponding electron gyro- frequency. In Eq. (10), the 
diffusion operator on the right hand side is raised to 2n. Here 
n is an integer and can take n = 1, 2, 3, • • •. The case n = 1 
stands for normal diffusion, while n = 2, 3, • • • corresponds 
to hyper- and other higher order diffusion terms. Eq. (10) 



alongwith Eq. (7) form a complete set of three dimensional 
compressible whistler wave model. For simulation purposes, 
we use two dimensional (2D) model by ignoring the variation 
in the z=direction and transforming the magnetic field as 
follows 

B(x,2/,t) = 2 x Vip(x,y,t) + 4>{x,y,t)z. 

Such representation preserves V-B = in 2D. Here if> and 4> 
are respectively orthogonal and logitudinal flux functions of 
the perturbed magnetic field B. This representation trans- 
forms Eq. (10) into its parallel and perpendicular compo- 
nents (of the magnetic field) as follows Abdalla et al (2001), 

A(0-d 2 V0) = ^ + [0,n + d 2 V 2 </.]-[^V 2 ^ + [/3,n],(ll) 



^(i> - d 2 V</>) = -\4>, 1> - d 2 V 2 </>] - \J3, v 2 # 



dt 



= -0,0, 



A 2 V 2 W> + /3), 



(12) 

(13) 
(14) 



where A = d e cj ce /u; 

The linearization of Eq. (10) about a constant magnetic 
field B = Boz + B, where Bo and B are respectively constant 
and wave magnetic fields, yields the following equation, 



w fc (l + d 2 fc 2 )B + -r^-ik k - B 



4.ne °- ^ 

On eliminating the wave perturbed magnetic field from the 
above relation, one obtains the following dispersion relation, 

dpk\\k 

(16) 



1+dlfc 2 ' 



k x + k y and k\\ = k-Bo. The use 



where Lo ce = eBo/m e c,k 2 
of Eq. (16) in Eq. (6) leads to the following relation between 
the wave magnetic field and the velocity field, 



±-k x B 

k 



(17) 



The rhs of Eq. (17), in combination with Eq. (6), cor- 
responds essentially to the whister wave perturbed veloc- 
ity field. This equation indicates that whistler waves consist 
of transverse fluctuations in the magnetized space plasma 
and they are produced essentially by rotational magnetic 
field that leads essentially to the velocity field fluctuations. 
On replacing the rhs in Eq. (17) with the perturbed ve- 
locity field, it can be shown that the whistler modes obey 
cquipartition between the magnetic and velocity field com- 
ponents as fc 2 |B| 2 ~ \V e \ 2 . The whistler wave activity can 
thus be quantified by how closely the characteristic modes 
obey the turbulent equipartition relation. In 2 and 3D cases, 
we have estimated this relationship respectively in Shaikh & 
Zank (2005) and Shaikh (2009) in the incompressible limit. 
Particularly interesting is the 2D case where Eq. (17) ex- 
hibits a linear relation between ipk and 4>k as k\ipk\ = \4>k\ 
(Shaikh & Zank 2005). It is further evident from Eq. (16) 
that there exists an intrinsic length scale d e , the electron 
inertial skin depth, which divides the entire turbulent spec- 
trum into two regions; namely short scale (kd e > 1) and long 
scale (kd e < 1) regimes. In the regime kd e < 1, the linear 
frequency of whistlers is ujk ~ k y k and the waves are dis- 
persive. Conversely, dispersion is weak in the other regime 
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Mode structures are shown in the low beta whistler turbulence simulations at a time when turbulence is fully saturated and is in the 
steady state. The effect of the background magnetic field is evident from these structures that are elongated along the x direction. The 
simulation parameters are: Box size is L x X L y = 2tt X 2n, numerical resolution is N x X N y = 512 X 512, electron skin depth is 
d e = 0.015 — 1.0, magnitude of constant magnetic field is So = 0.5, dissipation fi = 10 — 4 , time step dt = 10~ 4 . 



kd e > 1 since lou ~ k y /k and hence the whistler wave pack- 
ets interact more like the eddies of hydrodynamical fluids. 



3 NONLINEAR WHISTLER TURBULENCE 

We develop a two dimensional nonlinear fluid code to nu- 
merically simulate Eqs. (11) to (14) that describe low plasma 
beta whistler turbulence. The spatial descritization employs 
a pseudospectral algorithm (Gottlieb et al, 1977; Shaikh & 
Zank, 2006, 2007) based on a Fourier harmonic expansion of 
the bases for physical variables (i.e. the magnetic field, veloc- 
ity), whereas the temporal integration uses a Runge Kutta 
(RK) 4th order method. The boundary conditions are pe- 
riodic along the x and y directions in the local rectangular 
region of the solar wind plasma. The turbulent fluctuations 
are initialized by using a uniform isotropic random spec- 
tral distribution of Fourier modes concentrated in a smaller 
band of lower wavenumbers. While spectral amplitudes of 
the fluctuations are random for each Fourier coefficient, it 
follows a certain initial spectral distribution proportional 



to k~ a , where a is an initial spectral index. The spectral 
distribution set up in this manner initializes random scale 
turbulent fluctuations. We note that a constant magnetic 
field is included along the z direction (i.e. Bo = Box) to ac- 
commodate the large scale (or the background solar wind) 
magnetic field. 

The evolution of whistler fluid fluctuations are governed 
by the nonlinear mode coupling interaction processes. In the 
presence of a constant background magnetic field, turbulent 
fluctuations not only couple nonlinearly with each other, but 
they also propagate along the direction of the background 
magnetic field as small scale whistler wave packets. The in- 
teraction of whistler waves with turbulent fluctuations com- 
plicates the dynamical evolution. Additionally, by virtue of 
nonlinear interactions the larger eddies transfer their en- 
ergy to smaller ones through a forward cascade. According 
to Kolmogorov (1941), the cascades of spectral energy oc- 
cur purely amongst the neighboring Fourier modes (i.e. lo- 
cal interaction) until the energy in the smallest turbulent 
eddies is finally dissipated gradually due to the finite dis- 
sipation. This leads to a damping of small scale motions. 



© 0000 RAS, MNRAS 000, 000-000 



Inhomogeneous whistler turbulence in space plasmas 5 




Figure 2. 

Whistler turbulence in the large scale kd e < 1 regime exhibits a 
Kolmogorov-like inertial range power spectrum close to fc -7 / 3 . 

By contrast, the large-scales and the inertial range turbu- 
lent fluctuations remain unaffected by direct dissipation of 
the smaller scales. Since there is no mechanism that drives 
turbulence at the larger scales in our model, the large-scale 
energy simply migrates towards the smaller scales by virtue 
of nonlinear cascades in the inertial range and is dissipated 
at the smallest turbulent length-scales. A snap shot of fluc- 
tuations in density, magnetic and pressure fields is shown 
in Fig. (1). Consistent with 2D turbulence, a dual cascade 
phenomenon is observed Kraichnan (1965). In this process, 
the perpendicular component (0) of the magnetic field cas- 
cades predominantly towards the smaller scales whereas the 
parallel component (ip) exhibits large scales in its spectrum 
(see Fig la & lb). By contrast, density and pressure fields 
comprise of smaller scales (Fig lc & Id) as they follow a 
forward cascade. Since the spectrum of vorticity fields (V 2 tp 
and V 2 (f>) is dominated by the large k modes, they contain 
smaller scales (Fig le & If). 



4 ENERGY SPECTRA IN WHISTLER 
TURBULENCE 

The spectral transfer of turbulent energy in the iner- 
tial range is determined by neighboring Fourier modes in 
whistler turbulence. We find from our simulations that 
the mode coupling interaction follows a Kolmogorov phe- 
nomenology (Kolmogorov, 1941; Iroshnikov, 1963; Kraich- 
nan, 1965) that leads to Kolmogorov-like energy spectra. 
It is evident from Fig. (2) that whistler turbulence in the 
kd e < 1 regime exhibits a Kolmogorov-like k~ 7 ^ 3 spectrum. 
This inertial range turbulent spectrum, in the context of low 
beta whistlers, is further consistent with previous 2D work 
(Biskamp et al, 1996; Dastgeer et al, 2000a,b). Surprisingly, 
the density fluctuations do not modify the energy spectrum 
in the kd e < 1 regime. It turns out from the whistler wave 
dispersion relation that the wave effects dominate in the 
large scale, i.e. kd e < 1, regime where the inertial range tur- 
bulent spectrum depictes a Kolmogorov-like k~ 7 ^ 3 spectrum 
in our simulations. Our previous work, on the other hand, 



Figure 3. 

Spectrum of density fluctuations in whistler turbulence. The den- 
sity fluctuations are slaved to the magnetic field fluctuations and 
follow a Kolmogorov-like k~ 7 / 3 spectrum in the the kd e < 1 
regime. 

showed that the turbulent fluctuations in the smaller scale 
(kd e > 1) regime behave like non magnetic eddies of hydro- 
dynamic fluid and yield a k~ 5 ^ 3 spectrum (Shaikh 2009). 
The wave effect is weak, or negligibly small, in the latter. 
The observed whistler turbulence spectra in the kd e < 1 
regime in Fig. (2) can be followed from the Kolmogorov-like 
arguments (Kolmogorov, 1941; Iroshnikov, 1963; Kraichnan, 
1965) that describe the inertial range spectral cascades. We 
elaborate on these arguments to explain our simulation re- 
sults of Fig. (2) as follows. 

In the low plasma-/? regime, the whistler turbulence 
model described by the set of Eqs. (11) to (14) admits the 
following energy conservation law (Abdalla et al, 2001). 

E =\J dxdy [<j> 2 + (V^) 2 + dliV^f + dl(X7 2 4>) 2 

+d^{p 2 + V& 2 }]. (18) 

It is noted that the contribution due to the pressure and 
density fluctuations does not modify the energy conservation 
relation in the limit oj ce /cj pe < 1, B/Bo < 1 and /3 < 1. 

The kd e < 1 regime comprises the dispersive whistler 
waves and the total energy is dominated by first two terms 
in Eq. (18) such that E ~ <j> 2 + (Vt/>) 2 . Turbulent equipar- 
tition between the velocity and magnetic field fluctuations 
yields <f) ~ kip (Biskamp et al, 1996; Dastgeer et al, 2000a,b). 
Owing thus to the turbulent equipartition, the total en- 
ergy can be given as E ~ 4> 2 in the kd e < 1 regime. The 
group velocity of whistler waves in the kd e < 1 regime is 
v g ~ dui/dky ~ k ~ l~ . Assuming that the nonlinear 
transfer of energy in the inertial range is governed by the 
eddy interactions whose velocity is V e — z x V0. Apply- 
ing Kolmogorov-like dimensional arguments (Kolmogorov, 
1941), we obtain the k th Fourier component of the electron 
fluid velocity as ~ k(f>k. The convective time scales on 
which the eddies transfer energy in the inertial range can be 
estimated as r n i ~ l/(kvk) ~ l/(k 2 (f>k)- The nonlinear en- 
ergy cascade rates (e) are computed as e ~ Ek/r n i ~ k 2 4>\. 
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On using the Kolmogorov phenomenology that the 
spectral transfer is local and depends only on the energy 
dissipation rates and modes (Kolmogorov, 1941; Iroshnikov, 
1963), the energy spectrum can be given by E k ~ e c *fc /3 . 
Upon substituting the energy dissipation rates, we estimate 
the spectral energy as E k ~ e 2 ^ 3 k~ 7 ^ 3 . 

Thus, the energy spectrum k~ 7 ^ 3 derived on the basis 
of Kolmogorv-like arguments (Kolmogorov, 1941) is further 
consistent with our simulations in Fig. (2). 



5 EFFECT OF DENSITY FLUCTUATIONS 

We find from our simulations that the contribution due to 
the pressure and density fluctuations does not modify the 
energy conservation relation in whistler turbulence. Hence 
the density fluctuations do not influence the whistler wave 
cascades in the inertial range. This result appears to be a bit 
surprising at first, but our understanding is that the density 
fluctuations are only weakly coupled with the wave magnetic 
field and hence they are slaved to the magnetic field. During 
the evolution, density fluctuations simply follow the mag- 
netic field and exhibit a Kolmogorov-like energy spectrum. 
The density spectrum from our simulations is shown in Fig. 
(3) for the kd e < 1 regime. 

To quantitatively demonstrate that the density fluctu- 
ations couple only weakly with the whistler wave magnetic 
field, we develop a novel diagnostic to determine the contri- 
bution of density perturbation in the whistler waves. This 
originates essentially from the whistler wave relationship 
that relates the poloidal tj) and axial 4> components of the 
magnetic field flux functions and is described by the expres- 
sion k\ip k \ = \<f>k\ (Dastgeer et al, 2000a,b; Shaikh & Zank, 
2005). It should be further noted that this relation is valid 
for the incompressible whistler mode fluctuations. The effect 
of compressibility due to the density perturbations however 
enters through Eqs. (13) & (14). The whistler relationship 
is thus modified to include the density fluctuations and it 
reads as 



k 2 \tp k \ 2 



rik 



\ 2 k 2 



(19) 



where the left hand side corresponds to the energy asso- 
ciated with the wave field, whereas the right hand side of 
the expression describes the energy associated with the den- 
sity field. Based on the modified whistler relationship, we 
develop a parameter called as whistler parameter in the fol- 
lowing. 



xit) 



k 2 + 



(20) 



The whistler parameter is a dimensionless quantity that 
determines the contribution of density fluctuations in the 
whistler waves over the entire 2D turbulent spectrum. The 
physical picture emerging from this parameter can be de- 
scribed as follows; During the evolution of whistler waves 
and the corresponding density field, if there is a signicant 
amount of energy being transferred in the density fluctu- 
ations (i.e. comparable to the whistler wave energy), then 
the difference in the magnitude of the wave energy 
and density fluctuations \rik/X 2 k 2 — /3fc| 2 will be minimal. 
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Figure 4. 

Evolution of whistler parameter \ that determines the contribu- 
tion of density fluctuations in the whistler waves in the kd e < 1 
regime. The parameter \ approaching unity corresponds to the 
whistler wave dominated state in which density fluctuations are 
weakly coupled and they do not influence the energy cascade pro- 



Therefore a normalized contribution, obtained by dividing 
the difference by the total energy in the wave and the density 
field i.e. £ fcx *£ ky (k 2 \ipk\ 2 + \n k /X 2 k 2 - f3 k \ 2 ), will be much 
smaller than unity. This essentially corresponds to a state 
where x<Cl Such a criterion further characterizes a state in 
which the whistler wave magnetic field is greatly influenced 
or contaminated by the density fluctuations. By contrast, 
the limit x ~ ► 1 corresponds to a state in which density fluc- 
tuations contribute only weakly or negligibly small in the 
wave magnetic field. 

We follow the evolution of \ m our simulations. The 
result is plotted in Fig. (4). It is evident from Fig. (4) that 
X is increasing progressively and it is approaching unity i.e. 
X — ► 1. This means that energy associated with the wave 
field dominates over the energy in the density fluctuations. 
A decreasing trend of x> on l ne other hand, would have 
corresponded to state in which the energy in the density 
field was dominated over the wave. Since our simulations 
indiate that x ~ * 1, we believe that the effect of density 
perturbations is rather weak on the evolution of whistler 
waves and hence energy cascade rates are not altered. It 
is because of this reason that we find a Kolmogorov-like 
fc~ 7 / 3 energy spectrum in our simulations, a result similar to 
the one obtained previously by Biskamp et al (1996, 1999); 
Shaikh & Zank (2005); Cho & Lazarian (2004). 

Although Kolmogorov-like k~ 7 ^ 3 energy spectrum has 
been reported previously (Biskamp et al, 1996, 1999; Shaikh 
& Zank, 2005; Cho & Lazarian, 2004), they ignored the ef- 
fect of compressibility due to the density field in whistler 
turbulence. Our work described here is therefore different 
in a sense that we have included the density fluctuations in 
whistler turbulence and found that they do not influence the 
energy cascade processes in the inertial range spectrum. The 
inertial range turbulent spectrum in the kd e < 1 regime thus 
exhibits a Kolmogorov-like fc~ 7//3 power law. It thus follows 
from Figs (3) & (4) that the density field is simply slaved to 
the magnetic field and leads to Kolmogorov-like spectrum in 
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the wave-dominated kd e < 1 regime. Moreover, it should be 
noted from Fig. (4) that the evolution of whistler parameter 
is dominated entirely by the wave energy (over the density 
field) right from the outset. A primary reason of such domi- 
nance is ascribed to the form of Eq. (14) that relates density 
to the small scale density and potential fluctuations. This es- 
sentially means that the density field is dominated by the 
characteristic small scale fluctuations. Such small scale fluc- 
tuations tend to possess a weaker tendency to influence the 
turbulent spectrum in the kd e < 1 regime where the large 
scale whistler waves govern the entire nonlinear physics. This 
description further leads to another plausible explanation 
with regard to the time scales. The large scale energy con- 
taining whistler modes evolve on slower time scales [c.f. Eq. 
(16)]. By contrast, the small scale density fluctuations evolve 
on faster time scales. Owing to this temporal disparity asso- 
ciated with the two, the density field does not spend enough 
time with the wave field to influence its dynamical evolu- 
tion. Hence the former is substantially incapable of coupling 
with the wave field. This leads to the negligibly small con- 
tribution of the density fluctuations in the wave field. The 
inertial range cascade is therefore governed predominantly 
by the whistler mode interactions. 



6 ANISOTROPIC COMPRESSIBLE 
WHISTLER TURBULENCE 

We next quantify the degree of anisotropy mediated by the 
presence of large scale magnetic in the nonlinear 2D whistler 
turbulence. In 2D turbulence, the anisotropy in the k x — k y 
plane is associated with the preferential transfer of spectral 
energy that empowers either of the k x and k y modes. The 
background magnetic field is considered along the x direc- 
tion in our simulations. We therefore expect asymmetry in 
the evolution of the k x and k y modes. The anisotropy in the 
initial isotropic turbulent spectrum is triggered essentially 
by the background large scale magnetic field that regulate 
turbulent fluctuations to nonlinearly migrate the spectral 
energy in a particular direction. To measure the degree of 
anisotropic cascades, we employ the following diagnostics to 
monitor the evolution of feii mode in time. The fcii mode is 
determined by averaging over the entire turbulent spectrum 
that is weighted by k x which is aligned in the direction of 
the background magnetic field. 



(a) <j) field 



(b) v field 



*»(*) = 



'£J fc ttQ(M)l 2 
£JQ(M)l a 



(21) 



Here Q represents any of <f>, ip, V ip and V <j>. Similarly, 
the evolution of k± mode across the is background magnetic 
field determined by the following relation. 



kx(t) = 



£* 



\k y Q(k,t)\i 



£JQ(M)I 2 



(22) 



It is clear from these expressions that the k x and k y modes 
exhibit isotropy when k x ~ k y . Any deviation from this 
equality leads to a spectral anisotropy. We follow the evolu- 
tion of fey and kx modes in our simulations for long enough 
time. Our simulation results describing the evolution of fey 
and k± modes are shown in Fig. (5). It is evident from Fig. 
(5) that the initial isotropic modes k x ~ k y gradually evolve 




4 




time 



5 10 
time 



Figure 5. 



Evolution of anisotropic modes kn and kj_ averaged over the entire 
turbulent spectrum in low beta whistler turbulence. Initially, kn ~ 
fej_. Progressive development of anisotropy fcii < k± is ascribed 
to the presence of background magnetic field. 

towards an highly anisotropic state in that spectral transfer 
preferentially occurs in the k± mode, while the same is sup- 
pressed in fc|| mode. Consequently, the spectral transfer in 
k± mode dominates the evolution and the mode structures 
show elongated structures along the x-direction, see Fig. (1). 

The evolution of k± and kn shows a significant dis- 
parity in the two modes by virtue of an external magnetic 
field. The evolution of k± and fey is clearly different as the 
spectral cascade in the parallel wavenumbers is dramatically 
suppressed. The suppression is caused by the excitation of 
whistler waves, which act to weaken spectral transfer along 
the direction of propagation. This can be understood as fol- 
lows; We assume that the spectral transfer, essentially me- 
diated by propagating whistler waves in wavenumber space, 
can be described by a three wave interaction mechanism, for 
which the frequency and wavenumber resonance criteria are, 
respectively, expressed by 

and 

k 3 = ki + k 2 . 

The resonance conditions indicate that two whistler waves 
(o;i,ki) and (a^jka) mutually interact and give rise to the 
third wave (cj3, ka). Such conditions could, in principle, hold 
for a set of infinite waves as the indices T' and '2' are merely 
dummy indices. With the help of dispersion relation (say, 
in the kd e < 1 regime) and using the wavenumber k 3y — 



kly 

k-2„ 



+ k2 y , we can obtain 

_ k 2 + k 3 
ki — kz 



Let us now suppose the Kolmogorov turbulence hypothesis 
holds for EMHD nonlinear interactions as well, viz, that 
spectral transfer in the wavenumber space is local and occurs 
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efficiently amongst the most adjacent Fourier modes, i.e. for 
which \k — ki \ « |fci|. It, then, implies 

to, < ki v (23) 

since fci ~ k,2 « k$, thereby indicating that there is a very 
little cascading along the y-direction i.e. the magnetic field 
direction. Thus, the parallel wavenumbers (fey) appear to 
be suppressed and the spectral cascade mainly occurs in the 
perpendicular wavenumbers (k±). This, we suggest, explains 
the wavenumber disparity ({k±} =^ (fcy)) observed in our 
simulations [see Fig. (5)]. 



7 DISCUSSION AND CONCLUSIONS 

In summary, we present the results of freely decaying 
whistler turbulence calculated from an EMHD computer 
model. A major emphasis is to understand the effect of 
density fluctuations on the inertial range cascades in the 
whistler wave dominated kd e < 1 regime. In this regime, 
large scale whistler waves govern the energy cascade pro- 
cesses. Interestingly we find that despite strong density per- 
turbations, their effect on the cascade dynamics is incon- 
sequential. Hence they couple weakly with the large scale 
waves in the kd e < 1 regime. We find from our simulations 
that the density fluctuations do not influence the inertial 
range turbulent spectra. Consequently, the turbulent fluctu- 
ations in the inertial range are described by Kolmogorov-like 
phenomenology. Thus consistent with the Kolmogorov-like 
dimensional argument, we find that turbulent spectra in the 
kd e < 1 regime is described by k~ 7 / 3 . Our results are impor- 
tant particularly in understanding turbulent cascade corre- 
sponding to the high frequency (w > w c ») solar wind, space 
and astrophysical plasmas where characteristic fluctuations 
are comparable to the electron inertial skin depth. Our work 
might also be useful in describing why many space plasmas 
are described by by the Kolmogorov-like energy spectra de- 
spite the presence of strong density fluctuations and/or com- 
pressible effects. We find that density fluctuations exist on 
smaller scales and have typically higher frequency associ- 
ated with their evolution. Owing to these disparate length 
and time scales, they do not modify the large scale wave- 
dominated turbulent processes. It is this large scale dynam- 
ics in the kd e < 1 regime that leads to the Kolmogorov-like 
fc~ 7 / 3 spectrum in our simulations. 

Note that our simulations of inhomogeneous whistler 
turbulence pertains to decaying turbulence only. In princi- 
ple, turbulence can be driven. It remains to be seen whether 
the energy cascade processes in the driven whistler turbu- 
lence are altered by the density fluctuations. 
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